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DESIGN  STUDY  OP  DAMPING  TECHNIQUES 

I.  SCOPE  OP  STUDY 

Meetings  were  held  with  representatives  of  BuShips  Code  375  on 
30  June  and  27  August  at  which  the  scope  of  this  project  was 
discussed.  It  was  agreed  that  the  studies  under  the  present 
contract  should  cover  the  following  three  phases: 

Phase  I  -  Damping  by  Single  Homogeneous  Layers  and  by 
Simple  Constrained  Layers. 

A.  Analytical  studies  of  damping  on  large  plates, 
neglecting  the  effects  of  boundaries. 

B.  Analytical  studies  of  the  limitations  of  the 
results  of  Part  A,  including: 

1.  effects  of  boundaries 

2.  li.Mits  on  thickness  of  application 

3.  limits  of  temperature  and  frequency 

4.  effects  of  non-linearities 

5.  effects  of  water  load 

6.  effects  of  plate  curvature 

C.  Experimental  Studies 

Work  on  Part  A  of  this  phase  is  described  in  the  present  report. 
Work1  on  Parts  B  and  C  is  only  now  getting  underway. 
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Phase  II  -  Extensions  of  Analysis  to  Fibrous  Damping  Layers  and 
t-o  Stiff  Damping  Layers. 

A.  An  exploratory  analysis  of  the  constrained  fibrous 
damping  layer  aimed  at  understanding  the  mechanism 
of  damping  by  constrained  fibrous  damping  layers. 

B.  Analytical  study  of  constrained-layer  damping  when 
the  damping  layer  is  relatively  stiff.* 

It  is  planned  to  initiate  work  on  both  parts  of  this  phase  during 
the  October -December  quarter. 

Phase  III  -  Study  of  Combination  Damping  Treatments 

A.  Two  homogeneous  damping  layers 

or  B.  A  constrained  damping  layer  with  a  homogeneous 
layer  applied  over  it. 

It  was  agreed  that  the  decision  as  to  which  parts  of  Phase  III 
would  be  undertaken  would  be  postponed  until  work  in  Phase  II  was 
well  underway. 

During  this  quarter,  a  memorandum  was  submitted  by  the  contractor 
’outlining  work  that  has  been  and  is  being  undertaken  in  the  field 
of  vibration  damping  for  a  number  of  sponsors.  It  was  indicated 
that  although  these  studies  each  cover  different  damping  treat¬ 
ments  and/or  are  intended  for  different  applications,  each  should 
benefit  from  knowledge  gained  through  the  others. 
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II.  REPORT  ON  PHASE  I-A:  ANALYTICAL  STUDIES  OF  THE  DAMPING  OF 
LARGE  PLATES  BY  SIMPLE  VISCO-ELASTIC  DAMPING  TREATMENTS 

Introduction 

There  are  currently  two  distinct  methods  of  damping  the  flexural 
vibrations  of  plates  with  visco-elastic  materials.  The  older  and 
more  common  method  involves  the  application  of  a  relatively  thick 
layer  of  damping  material.  An  example  of  this  is  automobile 
undercoat.  The  equations  governing  this  type  of  damping  were  de¬ 
rived  independently  by  Oberst  and  Li^nardi^/.  The  damping 
mechanism  involves  the  internal  losses  of  the  visco-elastic  ma¬ 
terial  undergoing  alternating  extensional  strains.  The  damping 
achieved  is  roughly  proportional  to  the  product  of  the  resistive 
component  of  the  Young's  modulus*  and  the  square  of  the  thickness 
of  the  damping  layer. 

The  second  type  of  application  for  the  damping  of  plate  vibrations, 
called  a  "damping  tape"  or  M constrained  damping  layer”,  has  been 
developed  recently  in  the  aircraft  industry.  It  consists  of  a 
relatively  thin  layer  of  visco-elastic  material  and  a  stiff  cover¬ 
ing  ’'foil”.  The  damping  achieved  by  such  constrained  damping 
layers  involves  the  shear  motion  of  the  visco-elastic  material  be¬ 
tween  the  two  stiff  media.  The  damping  achievable  increases  with 
increased  stiffness  of  the  constraining  layer.  There  is  an  optimum 


*  This  component  is  sometimes  called  the  ’’loss  modulus”  to  dis¬ 
tinguish  it  from  the  elastic  component  or  "storage  modulus". 
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thickness  of  the  damping  material,  which  may  be  quite  small,  and 
which  is  a  function  of  the  temperature  and  frequency  range  for 
which  maximum  damping  is  desired.  The  equations  governing  the 
damping  by  constrained  damping  layers  were  recently  derived  by 
Kerwin^J^ 


A  more  general  analysis,  in  which  the  constrained -damping- layer 
configuration  is  treated  as  a  special  case,  has  been  carried 
out  by  BBN  and  presented  as  a  report  to  the  Office  of  Naval 
Research^/.  In  the  past  quarter,  we  have  taken  the  resultant 
equations  and  have  developed  charts  whereby  the  loss  factor  may 
be  calculated  from  knowledge  of  the  geometry  of  the  configuration 
and  of  the  elastic  moduli  of  the  materials.  We  have  then  pre¬ 
sented  these  results  in  the  form  of  charts  useful  in  understanding 
damping  results  and  in  designing  damping  treatments. 


Significance  of  the  Analyses 

The  analyses  of  the  damping  of  flexural  vibrations  by  visco¬ 
elastic  damping  layers  given  in  this  Progress  Report  are  ex¬ 
tensions  of  the  analyses  reported  in  Ref  7.  As  such,  they  are 
limited  to-  large,  flat  plates  for  which  the  effects  of  boundaries 
can  be  neglected.  None  of  the  secondary  effects  that  might  in¬ 
fluence  the  actual  performance  of  damping  treatments  on  submarine 
hulls  have  yet  been  included.  Specifically.,  the  effect  of  plate 
curvature,  external  pressure  and  of  water  load  have  not  yet  been 
considered.  Nor  have  the  practical  limitations  of  material 
properties  been  taken  into  account. 
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Despite  these  limitations,  we  believe  that  the  analyses  have 
appreciable  practical  value.  The  results  are  indicative  of  the 
maximum  damping  that  can  be  achieved  by  these  mechanisms.  The 
calculations  are  also  useful  in  comparing  various  damping  mecha¬ 
nisms  and  in  scaling  from  one  structure  to  another.  It  is  to  be 
expected  that  many  of  the  secondary  effects,  such  as  water  load 
and  plate  curvature,  will  have  similar  effects  on  different 
damping  treatments.  Thus,  many  comparisons  made  for  flat  plates 
will  also  apply  qualitatively  to  submarine  structures. 

No  experimental  confirmations  of  the  analyses  are  given  In  the 
present  report.  We  feel  that  the  experimental  confirmations 
already  reported  in  Refs  5  and  7  have  established  confidence  in 
the  basic  theories  involved. 


Homogeneous  Damping  Layers 

The  equations  governing  the  damping  of  homogeneous  visco-elastic 
applications  are  well  known,  having  been  drived  and  rederived 
many  times  (See,  for  example.  Refs  1,  3  and  7).  An  approximate 
relation,  which  holds  over  much  of  the  useful  range  is: 


where : 


^  =  14 


t)2e2  /  n2y 

Ei  [HiJ 


(1) 


tj  is  the  loss  factor* 

E  is  Young's  modulus 
H  is  thickness 

subscript  1  refers  to  the  base  plate 
subscript  2  refers  to  the  added  damping  layer 


*  The  relationships  of  tj  to  other  measures  of  damping  are 
in  Appendix  B. 
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This  equation  can  be  expressed  in  terms  of  the  weights,  W,  per 
unit  area  of  the  material  through  the  re^  >tion: 

W  =  ?H  (2) 

where  y  is  the  weight  per  unit  volume. 

The  resultant  expression  is: 


There  is  a  wide  variety  of  homogeneous -layer  damping  materials, 
having  best  damping  at  different  temperatures  or  having  different 
mechanical  properties,  e.g,,  adhering,  drying,  etc.  We  may  esti¬ 
mate  the  damping  achievable  with  the  better  products  when  applied 
to  steel  plate  by  taldtog  the  best  combination  of  parameters  that 
we  have  found  for  available  materials^/.  These  are: 

tj2e2  =  ^  x  dynes/cm2 
Y2  =  38.5  lbs/ft3 

=*  2  x  1012  dynes/cm2 
y,  =  480  lbs/ft3 
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The  resultant  expression  for  the  maximum  damping  factor  is: 


T)  ==  4 

'max 


W2N2 


Although  the  properties  of  the  various  compounds  vary  widely,  the 
greral  tendency  is  to  have  a  higher  resistive  stiffness  (loss 
modulus)  rE  associated  with  denser  materials,  so  that  the  coef¬ 
ficient  in  Eq  (4)  will  tend  to  remain  constant*.  Using  Eq  (4), 
we  have  constructed  Pig  1  to  show  how  the  approximate  maximum 
damping  varies  with  the  weight  of  the  application  relative  to 
that  of  the  plate.  In  Pig  2,  the  maximum  damping  Is  given  as  a 
function  of  the  weight  of  the  treatment  for  the  range  of  plate 
tVokriessrs  likely  to  be  encountered  In  submarine  construction. 

The  important  conclusion  from  the  foregoing  analysis  Is  that, 
although  the  relative  weight  of  the  damping  material  is  not  a 
fundamental  variable  in  the  damping  equations,  for  practical 
purposes  it  is  the  controlling  factor  In  determining  the  amount 
of  damping  that  can  be  achieved  with  homogeneous -layer  damping 
treatments . 

Constrained  Damping  Layers 

When  a  thin  layer  of  damping  material  is  constrained  by  a 
stiffener  layer,  the  resultant  damping  is  almost  entirely  attribu¬ 
table  to  the  shear  motion  of  the  damping  layer.  Analyses  of  this 


*  If  anything,  the  coefficient  will  be  somewhat  lower  with  denser 
materials. 
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mechanism  are  to  be  found  in  References  5  and  7.  Because  this 
theory  may  not  be  familiar  to  some  of  the  readers  of  this  report, 
a  brief  version  is  presented  in  Appendix  A.  The  final  theoretical 
equation  of  Reference  7  is  not  in  a  form  convenient  for  calculations; 
but,  it  is  possible  to  rearrange  the  terms  and  obtain  an  equation 
more  suitable  for  numerical  work.  This  has  been  done  in  Appendix  A. 
The  resultant  equation  is; 


Apr 


i  +  cr  +  d(i+ps)? 


(5) 


where;  p  is  the  loss  factor  of  the  shear  modulus  of  the 
damping  material.  (G*  =  Gg(l-fjp)). 


T  is  the  basic  parameter  governing  the  shear  damping. 

yon 

(r  =  -~r-  g-3-  for  steel  or  aluminum  plates.) 


.  A,  C  and  D  are  functions  of  relative  thicknesses  of 
the  foil,  damping  layer  and  base  plate  and  of  the 
materials  of  the  foil  and  base.  They  are  inde¬ 
pendent  of  the  properties  of  the  damping  layer. 


The  Appendix  also  contains  charts  whereby  A,  C,  D  and  r  may  be  cal¬ 
culated  from  the  material  properties  and  dimensions  of  the  layers. 
Other  charts  then  allow  the  calculation  of  the  resultant  damping 
factor.  Finally,  a  sample  calculation  is  carried  through  to 
illustrate  the  use  of  the  charts. 
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An  important  result  of  the  analysis  of  the  shear  damping  mecha¬ 
nism  is  the  dependence  on  the  properties  of  the  damping  layer  as 
given  by  the  shear  parameter,  P.  For  low  values  of  this  parameter, 
corresponding  to  a  damping  layer  with  little  resistance  to  shear 
motion.,  the  damping  is  relatively  low.  As  the  layer  becomes  stiffer, 
and  r  increases,  the  damping  increases.  The  damping  becomes  a  maxi¬ 
mum  at  a  moderate  value  of  r  (say,  0.01  <  r  <  0.1),  and  then  de¬ 
creases  with  further  increase  of  P,  as  the  layer  becomes  "too  stiff". 
The  peak  of  the  curve  is  relatively  broad,  the  loss  factor  dropping 
to  one-half  its  maximum  value  for  r  about  two  octaves  either  side 
of  the  optimum. 

The  shear  parameter  is  a  complicated  function  of  frequency  and 
temperature,  depending  as  it  does  on  the  shear  modulus  of  the  damp¬ 
ing  material,  which  is  a  function  of  both  frequency  and  temperature.  * 
In  addition,  r  depends  on  the  thickness  of  the  damping  layer.  For 
given  operating  temperatures  and  frequencies  it  is  usually  possible 
to  select  a  damping  material  and/or  layer  thickness  that  will  yield 
close- to-maximum  damping.  When  this  is  done,  the  damping  achieved 
is  primarily  a  function  of  the  stiffness  of  the  stiffener  layer 
relative  to  that  of  the  base  plate,  and  of  the  loss  factor,  p,  of 
the  damping-  material. 

To  estimate  the  damping  achievable  on  flat  plates  by  constrained- 
layer  damping  treatments,  we  have  assumed  that  the  proper  choice 
of  damping  material  has  been  made  to  give  optimum  damping  with  a 
layer  thickness  of  about  0.030  inches.  (This  layer  thickness  was 
chosen  as  the  minimum  that  can  safely  be  applied  to  heavy  steel 
plates  without  requiring  better'  than  normal  surface  finishes.) 

We  have  also  assumed  steel  as  the  stiffener  material,  although 
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aluminum  would  be  slightly  more  effective  on  an  equal  weight 
basis.  Assuming  values  of  B  of  0.5,  1  and  1.5,  we  have  plotted  , 
the  damping  factor  as  a  function  of  the  weight  of  the  treatment 
and  the  thickness  of  the  stiffener  layer.  The  results  are  pre¬ 
sented  in  Pigs  3-6  for  four  typical  thicknesses  of  steel  plate: 

1/4,  1/2,  1  and  2  inches.  The  three  curves  for  the  three  values 
of  B  may  be  interpreted  as  "readily  attainable",  "probable"  and 
"possible"  respectively. 

Comparison  of  the  Two  Types  of  Damping 

The  damping  achievable  by  homogeneous  damping  layers  depends 
roughly  on  the  square  of  the  weight  of  the  material,  while  that 
for  constrained  layers  depends  approximately  linearly  on  weight. 
When  we  compare  the  curves  of  Pig  2  with  the  corresponding  curves 
of  Pigs  3-6,  we  may  conclude  that  the  two  types  are  likely  to 
give  equal  values  of  maximum  damping  for  weights  between  10  and 
20  percent  that  of  the  base  plate. 

Below  10  percent,  the  constrained-layer  damping  will  probably  be 
the  more  effective;  while,  in  cases  where  more  than  20  percent  of 
the  plate  weight  can  be  devoted  to  damping,  then  the  homogeneous 
application  should  be  more  effective. 

An  important  corollary  of  this  result  is  that  different  types  of 
damping  treatments  cannot  be  directly  rank  ordered  if  they  are 
tested  at  ratios  of  treatment  to  plate  weights  different  from  those 
which  will  ultimately  be  used.  An  understanding  of  each  damping 
mechanism  is  necessary  before  meaningful  comparisons  can  be  made. 
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FIG.  I  MAXIMUM  DAMPING  OF  HOMOGENEOUS  LAYERS 
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FIG.  3  DAMPING  OF  1/4- INCH  STEEL  PLATES 

(10- POUND  PLATE) 
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FIG.  4  DAMPING  OF  1/2- INCH  STEEL  PLATES 

(20- POUND  PLATE) 
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THICKNESS  OF  STIFFENER  LAYER  IN  INCHES 

FIG.  5  DAMPING  OF  I -INCH  STEEL  PLATES 

(40- POUND  PLATE) 
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APPENDIX  A 

THE  DAMPING  OF  FLEXURAL  VIBRATIONS  BY  CONSTRAINED  DAMPING  LAYERS* 


Basic  Assumptions 

The  equations  that  are  presented  in  this  report  concern  the  damping 
attributable  to  the  shear  motion  of  a  thin  visco-elastic  layer  be¬ 
tween  two  stiff  plates.  Other  motions  could  occur  and  produce 
damping,  but  these  are  considered  negligible  relative  to  the  shear 
motion.  In  practice,  this  means  that  the  thicknesses  of  the  plates 
and  of  the  damping  layer  are  all  small  compared  to  the  shortest 
wavelengths  of  any  type  of  disturbance.  As  all  types  of  wave 
motion  are  ignored,  no  damping  factor  is  associated  with  the  mass 
per  unit  length  and  the  damping  factor  for  the  composite  plate  is 
simply  that  of  the  bending  rigidity.  All  effects  of  plate  curvature 
and  of  water  load  are  ignored  in  this  elementary  analysis. 

Flexural  Rigidity  of  Composite  Plate 

Figure  A-l  represents  an  element  of  a  three-layer  structure  in 
flexure,  showing  the  shear  displacement  of  the  middle  layer. 

Layer  1  represents  the  original  undamped  plate,  the  second  layer 
is  the  damping  material,  and  the  third  layer  the  stiff  constraining 
layer.  The  angle  <t>  is  the  flexural  angle  of  the  element  and  the 
shear  strain,  f,  is  measured  relative  to  $  in  the  opposite  direction. 
The  x-dlrection  is  chosen  as  the  direction  of  propagation  of  the 
straight-crested  flexural  wave.  It  is  assumed  that  all  of  the 
layers  vibrate  in  phase  and  that  the  flexural  angle  of  an  element 
of  foil  Is  the  same  as  that  for  the  corresponding  element  of  the 
base. 


♦The  equations  are  derived  in  more  detail  in  Ref  7.  Some  of  the 
figures  used  here  come  from  that  report. 

A-l 


20 


Report  No.  591 


Bolt  Beranek  and  Newman  Inc. 


Figure  2  shows  the  various  thicknesses  and  distances  used  in  the 
analysis.  The  displacement  of  the  neutral  plane  is  represented 
by  D. 

The  distribution  of  extensional  strain  and  stress  for  the  three- 
Jayer  element  is  illustrated  in  Fig  3*  It  is  assumed  that  the 
stiffness  of  the  damping  layer  is  small  compare^  tc  that  of  the 
other  two  layers,  so  that  the  net  extensional  stress  of  the  damping 
layer  can  be  neglected.  As  the  net  force  on  the  composite  element 
is  zero,  the  individual  layer  forces  of  the  top  and  bottom  layers 
are  equal  in  magnitude  and  opposite  in  direction.  Letting  the 
force  on  the  constraining  layer  be  F,  we  may  write; 

P  =  K3(H31-D)  *  -  *3^  !£  "  hD  H  (A-1' 


fch 

where  is  the  stiffness  of  a  unit  area  of  the  1  layer. 

The  bending  moment  about  the  neutral  plane  of  the  composite  struc¬ 
ture  is  the  sum  of  the  moments  of  the  individual  layers  about 
their  own  centers  and  of  the  products  of  the  force  on  each  layer 
by  the  distance  of  the  center  of  that  layer  from  the  neutral  plane 
of  the  composite.  The  expression  for  the  bending  rigidity  is 
thus : 


B  *  B$73x  =  15  K1H12  +  T2  k3H32  +  K1D  H31  (A"2 

3$ie  problem  now  reduces  to  finding  the  displacement  of  the  neutral 
plane  of  the  composite  relative  to  that  of  the  base.  Solving 
Eq  (A-l)  for  D:  : 
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K, 


D  = 


K1  +  K3 


-  H 


2  W 


It  is  seen  that  this  displacement  is  a  function  of  the  relative 
strain  of  the  damping  layer. 


The  shear  strain  of  the  middle  layer  is  related  to  the  shear  stress 
and  the  shear  modulus.  As  shown  in  Pig  4,  the  shear  force  is  pro¬ 
portional  to  the  space-derivative  of  the  force  experienced  by  the 
foils 


*  = 


-1  oF 
G2  kx 


(A-4) 


where  G2  is  the  shear  modulus  of  the  damping  layer.  Except  near 
the  edges  of  the  structure,  the  vibratory  motion  is  characterized 
by  expressions  of  the  forms 


+JPnX  Jo>nt 
e  e 


(A-5) 


where  p_  is  the  wave-number  of  the  n4'  mode  and  cd>  is  the  corres- 
n  n 

ponding  angular  frequency.  Ic  follows  that  the  second-derivative 
of  the  shear  strain  is  related  to  the  strain  itself  by: 


a2*  _ 


pn2* 


(A-5) 


A-3 


•?  o 

u  u 


and  it  follows  that 
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(aVa*2)  _ 

(a'V^) 


(A-7) 


The  quantity  In  parentheses  is  a  dimensionless  coefficient  that  ex¬ 
presses  the  properties  of  the  shear  layer  relative  to  the  stiffness 
of  the  basic  structure  and  the  wave  number  of  the  vibration.  This 
is  a  fundamental  parameter  of  the  damping  tape  problem.  Putting 


K3.H2P 

substituting  Eq  (A-7)  in  Eq  (A-3)>  solving  for  the  displacement  D 
and  then  using  Eq  (A-2),  the  final  expression  for  the  bending 
rigidity  is: 


B  -  ^  KlHl2 


+  15-  K3«; 


ralK3H3! 


k3  + 


TTK^ 


(A-9) 


The  dissipative  properties  of  the  middle  layer  are  considered  by 
expressing  the  shear  modulus  as  a  complex  quantity: 

G2*  =  G2(l  +  jp)  (A-10) 

where  both  G0  and  p  are  functions  of  temperature  and  frequency. 
Substituting  Eq  (A-10)  for  G*  into  Eq  (A-8)  for  r*  and  then  solving 
Eq  (A-9)  for  B*  one  finds: 


A-4 
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12  KlHl  +  12  K3H3 


The  loss  factor  for  the  composite  structure  is: 


(A-ll) 


pr  Kih31' 


T5k1H1+  T5k3h3  1+2r(”^+rS('RJ‘)  (1+e2)+nc1H|1  1+r(^~)f+@2') 


(A-12) 


Porm  for  Calculations 


Equation  (^12)  for  the  loss  fa -tor  of  a  plate  with  constrained  layer 
damping  can  be  rearranged  to  facilitate  calculations.  The  major 
independent  variable  is  the  shear  parameter  r,  while  the  relative 
tnicknesses  of  the  damping  layer  and  of  the  constraining  layer,  as 
well  as  the  stiffness  of  the  "foil”  and  the  loss  factor  of  the 
damping  material,  are  parameters.  First,  all  properties  of  the 
added  layers  can  be  normalized  by  dividing  by  the  corresponding 
property  of  the  base  plate: 


h3  =  H3^a 


(A-13) 


A  EA 

k3  ~  V*!  "  =  e3h3 


(A-H) 
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The  distance  between  the  centers  of  the  constraining  layer  and 
the  base  plate  is  normalized  by  dividing  by  the  thickness  of  the 
base  plates 


31 


«3l 

HJ” 


i  H,  +  H„  +  i  -i 

~  ~~  Hr~  ~  ~  - 1 <1+2W 


(A-15) 


This  distance  appears  several  places  in  the  calculations,  always 
in  combination  with  another  term;  so,  we  define: 


A  = 


12h~, 

1  +  k^h^ 


3(l+2h24h3)J 


(A-16) 


For  reasons  which  will  become  clear  later,  we  define  two  additional 
coefficients  that  are  dependent  only  on  hg,  h^  and  e^  (or  k^): 


C  5  A  +  2 


a  +  k. 


(A-17) 


rl  +  k 


3 


k- 


1  +  Ai 


1  +  k- 


(A-18) 


With  these  two  additional  abbreviations,  the  equation  for  the  loss 
factor  reduces  to: 


n  = 


Mr 


i  +  cr  +  d(i  +  p2)  rc 


(A-19) 


which  is  a  form  that  is  suitable  for  basic  calculations. 

A-6 
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Optimum  Damping 

The  loss  factor  for  a  plate  with  a  constrained  damping  layer 
treatment  is  highly  dependent  on  the  shear  parameter  of  the  damping 
layer.  At  low  values  of  r  it  is  proportional  to  1%  while  at  high 
values  it  varies  inversely  with  P.  There  is  a  value  of  r  that  gives 
maximum  damping.  If  we  call  this  value  rop^.,  we  may  solve  for  it 
by  taking  the  derivative  of  Eq^-19)  and  setting  this  equal  to  zero. 
The  result  is: 


opt 


+  p2 


(A-20) 


Equation  (A-19)  for  the  loss  factor  can  now  be  put  into  its  most  simple 
form  by  defining  three  reduced  coefficients: 


2>Jd 


(A-21) 


2  Jd 


(A-22) 


P'  =  tan"1  p 


(A-23) 


The  resultant  expression  for  q  is: 


R  = 


2A«smp'(r/ropt) 


1  +  2c1  cos  p«  {P/ropt)  +  (r/ropt) 


(A-24) 


A-7 
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where 


opt 


cos  0* 


(A-25) 


The  maximum  damping  achievable  with  a  given  damping  treatment  is 
given  by ; 


TWx 


a  A* 


sin  p* 

1  +  C'  cos  p1 


A‘f(p) 


(A-26) 


The  maximum  damping  is  a  function  of  the  loss  factor  of  the  damping 
material  as  well  as  of  the  relative  dimensions  of  the  two  layers  of 
the  treatment.  At  non-optimum  values  of  r,  the  damping  factor  is 
less  than 


W  (r/ront  “  15 

1  +  S(1  +  C  cos  (i')(r/ropt) 


(A-27) 


This  function  is  geometrically  symmetric  with  respect  to  r/ropf 


Graphs  of  General  Functions  Used  in  Damping  Calculations 

The  general  results  of  the  analysis  of  damping  by  constrained  damping 
layers- are  expressed  by  Eqs  (A-25), '  (A-r26),  and  (A-27)  in  terms  of  • 
coefficients  A',  C*  and  p’  which  are  functions  of  the  physical 
properties  of  the  layers  making  up  the  treatment.  The  important 
parameter  is  k^,  the  relative  stiffness  of  the  outer  layer.  Other 
parameters  are  hg  and  h^,  the  latter  being  related  to  k^  through 


A-8 
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=  e3^3*  To  enable  rapid  calculations  of  damping  and  the  later 
development  of  design  charts,  graphs  are  needed  for  the  various 
coefficients  in  terms  of  the  basic  physical  quantities. 

For  making  our  calculations,  we  have  chosen  values  of  e^  of  1/3, 

1  and  3.  If  both  the  foil  and  the  plate  are  made  of  the  same 
material,  then  e^  *  1.  If  one  is  of  a  light  metal  and  the  other 
a  heavy  metal,  then  the  ratio  of  the  Young’s  moduli  will  be  close 
to  1/3  or  3,  depending  on  which  one  is  on  top.  It  is  unlikely 
that  the  damping  layer  will  have  a  thickness  greater  than  a  third 
that  of  the  plate,  so  values  of  hg  =  0,  0.1,  0.2  and  0.3  have  been 
selected. 

Figures  5  to  7  give  the  dependence  of  ropt  on  e^,  hg  and  k^- 
Actually,  the  three  sets  of  curves  are  very  similar  for  k^  <  0.1, 
only  diverging  for  the  higher  values.  Within  the  accuracy  of  many 
calculations.  Fig  6  can  be  used  for  all  values  of  By  and  the  curve 
for  h2  »  0.1  is  representative  of  all  values  of  hg.  The  function 
of  cos  0'  needed  to  complete  the  calculation  of  r  is  given  in 
Fig  12  as  a  function  of  0. 

The  maximum  loss  factor  is  a  function  of  £  and  of  the  parameters 
e^*  hg  and  k^.  Figures  8  to  10  show  its  dependence  on  the  latter 
group.  In  this  case  variations  of  both  hg  and  e ^  are  important. 

The  function  f(0)  is  a  function  of  C'  as  well  as  of  0.  Ac  shown 
in  Fig  11 ,  C '  Is  a  3lowly  varying  function  of  A',  having  values 
between  1.0  and  1.4.  (That  C’  should  be  related  to  A’  was  shown 
both  analytically  and  empirically  before  Fig  11  was  drawn.)  We 
may  therefore  relate  all  functions  of  C'  to  functions  of  A', 
choosing  three  values  of  C  as  representative  of  the  full  range 
of  A'  as  follows: 

A -9 
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0  <  A‘  <  .3  C*  =  1.02 

.3  <  A‘  <  .6  C‘  »  1.10 

.6  <  A'  <  1  C'  «  1.26 

In  Pig  12,  the  various  functions  of  0'  are  plotted  against  0,  using 

the  three  representative  values  of  C*  where  required.  It  is  clear 
that  the  variations  caused  by  A';  i.e.,  by  k^,  e^  and  hg,  are  small. 
As  the  loss  factor  of  the  material  is  usually  not  known  with  any 
great  precision,  it  is  reasonable  to  neglect  the  small  ex’fect  of 
these  other  factors. 

Dependence  on  T/r^. 

The  ratio  of  the  damping  for  any  value  of  r  to  the  maximum  achievable 
at  rQpt  is  primarily  a  function  of  the  ratio  of  r  to  Topt  and  second¬ 
arily  a  function  of  0  and  the  dimensions  of  the  layers.  We  may 
neglect  the  variation  with  layer  dimensions  and  consider  only  the 
effect  of  0.  Even  this  is  a  slow  variation  and  only  a  few  typical 
values  of  (1  +  C1  cos  0’)  need  be  used  in  £q  (A-27).  The  values  used 
are  the  following: 

0  <  0  <  .6  (1  +  C'  cos  0')  *  2.00 

.6  <  0  <  1.4  »  1.75 

1.4  <  0  <  3  =  1.50 

Figure  13  shows  how  varies  with  r/Topt  for  these  three 

ranges  of  0.  Because  of  the  geometric  symmetry  of  the  function 
one  set  of  curves  applies  whether  the  shear  parameter  be  greater 
or  lesa  than  optimum.  As  zhe  variation  with  0  is  relatively 


A- 10 


29 


Report  No.  591 


Bolt  Beranek  and  Newman  Inc. 


small,  and  as  0  is  most  often  between  0.6  and  1.4  for  typical 
damping  materials,  it  will  usually  be  within  the  desired 
accuracy  simply  to  use  the  middle  curve.  To  within  +  !(#: 


_!L__ 

TWx 


3-65  (r/ropt) 

(i  +  rAopt)* 


(A-28) 


no  matter  what  the  value  of  the  loss  factor  0. 


The  Shear  parameter 

The  shear  parameter  defined  by  Eq  (A-8)  is  a  function  of  the  shear 
modulus  and  thickness  of  the  middle  layer,  of  the  stiffness  of  the 
base  plate,  and  of  the  wave-number  of  the  vibrations.  The  shear 
modulus  and  wave-number  are  both  functions  of  frequency.  The 

former  is  a  property  of  the  viscous  material,  while  the  wave-number 

-  »*■ 

is  related  to  the  frequency  through  the  velocity  of  the  flexural 
waves  by: 


The  flexural  velocity  is  influenced  by  the  added  stiffness  and 
mass  of  the  damping  layer.  If  we  assume  that  this  is  a  second- 
order  correction,  then  we  may  estimate  Cg  from  the  properties  of 
the  undamped  plate.  The  formula  for  the  flexural  wave  velocity 
for  an  undamped  plate  is: 
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where  co  is  the  usual  velocity  of  sound  in  the  plate.  Using 
this  relationship,  the  expression  for  the  shear  parameter 
becomes : 


<*£ 2f 


HT  1C 


(A-31) 


For  most  metal  plates,  the  longitudinal  velocity  of  sound  is 
closely  17,000  ft/sec,  so  that  the  shear  parameter  can  be 
calculated  from: 


,  790  °2  .  720  °2 


(A-32) 


where  is  the  Young's  modulus  of  the  plate  material  and  a 
correction  has  been  made  for  Poisson's  ratio. 


11  p 

The  value. of  Young's  modulus  is  about  7.3  x  10  dynes/cm  for 

11  p 

aluminum  and  20  x  10  dynes/cra  for  steel.  The  shear  modulus 

2 

of  visco-elastic  materials  varies  from  10  dynes/cm  when  very 

10  p 

soft  to  about  10  dynes/cra  when  hard.  The  thickness  of  the 
damping  layer  could  be  as  small  as  1C~^  ft  or  as  large  as 
2  x  10~2  ft.  It  follows  that  possible  values  of  T  could  vary  by 
as  much  as  10^,  at  a  given  frequency.  However,  significant 
damping  can  only  be  achieved  in  the  range:  0.002  <  r  <  2.  By 
assuming  values  of  Yeung's  modulus  appropriate  to  the  plate 
material,  we  can  plot  r  as  a  function  of  the  shear  modulus  of* 
the  damping  material,  with  frequency  and  damping  layer  thickness 
as  parameters.  Fig  14  is  such  a  plot  to  be  used  to  estimate  the 
shear  parameter  when  the  base  plate  is  made  of  aluminum,  while 
Fig  15  applies  to  steel  plates. 

A-12 
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Use  of  the  Charts 

The  charts  presented  in  Pigs  5  through  15  can  be  used  to  calculate 
the  shear  damping  of  a  given  constrained  damping- layer  treatment. 
Thus,  knowing  the  thicknesses  of  the  foil  and  damping  layer  rela¬ 
tive  to  that  of  the  base  plate,  one  can  readily  compute  h^  and  hg. 
The  materials  of  the  plate  and  foil  determine  e^  and  thus  k^. 

Going  into  the  appropriate  char.t  among  Pigs  5  through  7 ,  we  obtain 
ropt/cos  0' .  In  the  same  way.  Pigs  8  through  10  give  A‘  *= 
t^^/ffp).  These  two  results  are  independent  of  the  damping  ma¬ 
terial  and  therefore  of  temperature  and  frequency.  Given  the 
properties  of  the  damping  material,  the  appropriate  values  of  the 
loss  factor  can  be  tabulated  for  each  combination  of  frequency 
and  temperature.  Pig  12  gives  both  cos  0'  and  f(p)  for  each  value 
of  0.  One  can  then  calculate  rjmax  and  r  ^  for  each  combination 
of  frequency  and  temperature.  It  then  only  remains  to  determine 
the  actual  values  of  T  from  Pigs  14  or  15  and  to  calculate  the  re¬ 
duction  in  tj  caused  by  nop-optimum  values  of  r,  using  Pig  13. 

As  an  example  of  the  use  of  these  charts  let  us  carry  through  the 
calculation  of  the  damping  curve  for  a  damping  treatment  consisting 
of  1/32  inch  polyisobutylene  and  1/8  inch  3teel  on  a  1-inch  steel 
plate  at  room  temperature.  As  the  constraining  layer  is  the 
same  material  as  the  plate,  =  1.  From  geometry:  • 


>2  =32 
1 


=  .030 


*3  *  h3  =  "  *125 
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Prom  Pig .  6: 


-r°pfc  =  JL 

cos  p* 


.091 


Prcm  Pig  9: 


.195 


Prom  Pig  11:  c*  *=  1.02 


Prom  references  8  and  9  (or  Fig  14  of  Ref  7)  concerning  the  dynamic 
mechanical  properties  of  poly isobutylene,  we  can  determine  the 
shear  modulus  and  loss  factor  for  several  typical  frequencies  at 
room  temperature : 


Frequency 

jG 

B. 

30  cps 

4.2  x  10^  dynes/cm2 

0.6 

.  100  cps 

6.0 

0.85 

300  cps 

9.2 

1.2 

1000  cps 

16.5 

1.5 

From  Pig  12,  we  obtain  cos  and  f(p)  and  thus  calculate  rQpt  and 

w 


Frequency 

cos 

ropt 

gei 

TWx 

30  cps 

0.85 

.077 

0.28 

.055 

100  cps 

0.76 

.069 

0.37 

.072 

300  cps 

0.65 

.059 

0.47 

.092 

1000  cps 

0.55 

.050 

0.53 

.104 
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Prom  Pig  15 ,  the  shear  parameter  can  be  read  and  the  ratio  to 
optimum  calculated: 


Frequency 

fH2 

£2 

r 

JZSopt 

lopi/T 

30  cps 

0.9 

4.2  x  10 

.021 

0.27 

3.7 

100  cps 

-m* 

6.0 

.0085 

0.123 

8.1 

300  cps 

9 

9.2 

.0044 

0.075 

13.5 

1000  cps 

30 

16.5 

.0024 

0.048 

21 

Finally,  the  reduction  in  tj  caused  by  non-optimum  r  is  obtained 
from  Pig  13  and  tj  then  computed: 


Frequency 

_JL_ 

30  cps 

0.64 

.035 

100  cps 

0.36 

.026 

300  cps 

0.225 

.021 

1000  cps 

0.14 

.0145 

A- 15 
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FIG.  A- 1  ELEMENT  OF  A  THREE-LAYER  PLATE  IN 
FLEXURAL  VIBRATION,  SHOWING  THE 
FLEXURAL  ANGLE  <P  AND  SHEAR  ANGLE  f 
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FIG.A-2  DIMENSIONS  USED  IN  ANALYSIS  OF 
A  THREE- LAYER  PLATE  IN  FLEXURE. 
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FIG.A-3  EXTENSIONAL  STRAIN  AND  STRESS  DISTRIBUTIONS 
FOR  THREE -LAYER  PLATE  ELEMENT  IN  FLEXURE. 
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FIG.  A- !2  FUNCTIONS  OF  £ 
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FIG.A-15  SHEAR  PARAMETER  , ON  STEEL  PLATES 
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APPENDIX  B 

MEASURES  OF  DAMPING  EFFECTIVENESS 


Loss  Factor,  q 

The  velocity  of  propagation  of  flexural  waves  on  an  undamped  and 
unloaded  plate  is: 

4  — 

CB  =  ^  Jl  (B-1) 

where  B  is  the  flexural  rigidity  per  unit  width  and  m  is  the  mass 
per  unit  area.  As  the  velocity  is  dependent  on  the  frequency, 
flexural  waves  are  said  to  be  dispersive.  We  may  treat  the  damp¬ 
ing  caused  by  the  internal  friction  of  the  materials  composing  the 
plate  by  defining  a  complex  flexural  rigidj 

B*  =  B'  +  JB"  =  B(1  +  jTjg)  (B-2) 


where  is  the  loss  factor.  If  there  is  damping  caused  by  radi¬ 
ation  or  other  wave  effects,  then  the  mass  is  expressed  as  a  complex 
quantity: 


m*  =  m'  -  jm"  -  m(l  -  jr}^) 


(B-3) 


If  both  types  of  damping  occur,  then: 


(1  +  Jt)) 


(B-U) 


B-l 
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where  the  effective  ratio  of  the  rigidity  to  mass  is 


1  -  Vto 


/B\  _  /B\  _  'B 

We  W  1  +  nf; 


% 


(B-5) 


and; 


r,B  +  1m 

1  1  -  vw 


(B-6) 


As  the  individual  loss  factors  are  seldom  greater  than  0.2,  the 
total  effective  loss  for  the  plate  is  closely  equal  to  their  Siam. 
If  the  total  loss  factor  is  less  than  0.^,  th(  n  the  expression 
for  the  complex  velocity  of  propagation  is  i> accurately; 


°i 4  ^t/I  (i  +  i  3) 


(B-7) 


The  loss  factor,  tj,  defined  by  Eq  (B-4)  is  a  common  measure  of  the 
effectiveness  of  damping. 

Figure -of -Merit,  Q 

The  frequency  response  of  a  finite  plate  contains  many  resonances. 
Another  measure  of  damping  is  the  figure-of -merit  or  Q  of  the 
resonances,  defined  by; 


«s-s? 


(B-8) 


FJ-2 
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where  Af  Is  the  bandv?idth  between  the  two  half-power  points  cf 
the  resonance.  The  Q  of  a  resonance  is  equal  to  the  reciprocal 
of  the  loss  factor  at  the  frequency  of  the  resonance: 

Q  =  i  (B-9) 
provided  the  damping  factor  is  not  greater  than  about  0.2. 


Fraction  of  Critical  Damping 

In  many  vibratory  problems  in  mechanics  the  free  motion  is  governed 
by  a  linear,  second-order  differential  equation  of  the  form: 

mx  +  cx  +  kx  =  0  ^B-10) 


Only  if  the  damping  coefficient,  c,  is  less  than  a  critical  value, 

cr  =  2  \ZnflT  (B-ll) 

is  the  motion  oscillatory.  The  ratio  of  the  damping  coefficient, 
c,  to  the  critical  damping  coefficient  is  then  a  measure  of  the 
damping  effectiveness. 


For  low  values  of  damping  the  shapes  of  the  resonance  curves  of  a 

plate  resonance  and  a  lumped -parameter  system  whose  notion  is 

governed  by  Eq  (B-10)  are  similar.  Therefore  a  value  of  c/c.  can  ' 

c 

be  attributed  to  a  given  value  of  the  loss  factor  tj.  The  Q  of  a 
simple  mechanical  system  is  given  by: 


Q  = 


(B-12) 
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whence  the  equivalent  fractional  critical  damping  of  a  plate 
vibration  is: 


c 


(b-13) 


Equation  (B-13)  can  be  used  as  the  definition  of  the  fractional 
critical  damping  for  flexural  vibration*  of  plates. 


Space  and  Time  Decay  Rates 

Damping  may  b^  measured  by  observing  the  time -rate -of -decay  of  a 
vibration  after  the  driving  force  is  removed,  or  by  measuring  the 
spatial  decay  when  a  steady-state  force  is  applied  at  one  place. 

For  values  of  q'less  than  0.2,  these  decay  rates  can  be  expressed 
in  terms  of  rj  by  the  following: 

Dt  =  27.3Rf0  (db/sec)  (B-14) 

=  13 . 6q  (db/wave length)  (B-15) 

*s  13. 6rj/^  (db/unit  length)  (B-l6) 

Significance  of  the  Loss  Factor 

A.'.'.-ough  any  of  the  many  ways  of  describing  damping  may  be  used  for 
small  amounts  of  damping,  it  is  best  to  express  analytical  results 
in  terms  of  a  parameter  that  has  meaning  for  highly  damped  as  well 
as  lightly  damped  structures  and  which  applies  tr  distributed 
systems  as  well  as  lumped-parameter  systems.  The  loss  factor-  rj, 
defined  by  Eq  (B-4)  is  quite  general.  For  this  reason,  we  shall 
express  all  of  our  results  in  terms  of  q,  and  the  other  quantities 
may  then  be  estimated  from  Eqs  (B-9)>  (B-13)>  (B»l4),  (B-15)>  and 
(B-36). 


B-4 


